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ABSTRACT. Thermal ions in spherical tokamaks 
have two adiabatic invariants: the magnetic moment and 
the longitudinal invariant. For hot ions, variations in 
magnetic-field strength over a gyro period can become 
sufficiently large to cause breakdown of the adiabatic 
invariance. The magnetic moment is more sensitive 
to perturbations than the longitudinal invariant and 
there exists an intermediate regime, super-adiabaticity, 
where the longitudinal invariant remains adiabatic, 
but the magnetic moment does not. The motion of 
super-adiabatic ions remains integrable and confinement 
is thus preserved. However, above a threshold energy, 
the longitudinal invariant becomes non-adiabatic too, 
and confinement is lost as the motion becomes chaotic. 
We predict beam ions in present-day spherical tokamaks 
to be super-adiabatic but fusion alphas in proposed 
burning-plasma spherical tokamaks to be non-adiabatic. 

Consider a charged particle, with mass m and 
charge q, gyrating with speed v in a magnetic field, 
with field strength B. Assume axi-symmetric toroidal 
geometry with coordinates (R, </>, Z). The motion 
preserves the energy E (because the Lorentz force 
performs no work) and canonical toroidal momen- 
tum i-0 (because the toroidal coordinate is ignor- 
able). There is no third exact constant of motion, 
but the normalized magnetic moment 
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where v± is the velocity component perpendicular to 
the magnetic field (vj_ = v 2 — v 2 , V\\ = v-B/JE?) and Bo 
is the on-axis magnetic field, is an adiabatic invari- 
ant, i.e. as long as the variation in the magnetic field 



experienced by the particle during one gyro period 
is small, A oscillates at the gyro frequency around a 
constant average. If we introduce the adiabaticity pa- 
rameter, e, the condition of slow (adiabatic) variation 
becomes 

where the gyro radius g = vj_/Q, with the gyro 
frequency Q = \q\B/m. The approximate constancy 
of the magnetic moment was first pointed out 
by Alfven §. Just as the magnetic moment is 
associated with the gyro motion, a second adiabatic 
invariant is associated with the slower drift motion: 
the longitudinal invariant J|| = §v\\M 

When £ <C 1, there are thus three constants of 
motion, one for each degree of freedom, i.e. the mo- 
tion is integrable, the orbits are quasi-periodic and, 
in the absence of collisions, a particle will be eter- 
nally confined. As the perturbation grows and e < 1, 
a dramatic change in the character of the motion, 
from bounded to unbounded, will occur, as the adi- 
abaticity breaks down, the third constant of motion 
is lost, the motion becomes non-integrable, and the 
orbits become chaotic. This process was first inves- 
tigated numerically Q and later fully explained by 
the rigorous KAM (Kolmogorov, Arnold & Moser) 
theory Subsequently, Chirikov provided an ap- 
proximate, but more practical, approach when he in- 
troduced his semi-empirical resonance-overlap crite- 
rion and the standard mapping |). In the Chiriko- 
vian model, the breakdown of adiabaticity is caused 
by strong nonlinear resonance between the periodic 
motions occurring on different timescales. 
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The compact and low-field spherical tokamak was 
introduced by Peng and Strickler [0] . The first exper- 
iment, START §, showed encouraging results and 
has been succeeded by two larger machines, MAST || 
and NSTX |1. Larger, burning-plasma, experiments 
have also been proposed p0| . The B contours of 
spherical tokamaks are characterized by a magnetic 
well, i.e. a local minimum in B, centered outside of 
the magnetic axis [see Fig. [j]]. For the case shown 
here, B min = 0.21 T and B max = 1.26 T. The low 
value and large variation of B in the plasma (a factor 
of six) leads to a low threshold energy for the onset 
of non-adiabatic motion. Furthermore, the threshold 
energy will vary substantially over the cross section. 
We will demonstrate the consequences of this by solv- 
ing the Lorentz equation of motion for 80 kV deu- 
terium ions moving in the NSTX magnetic equilib- 
rium shown in Fig. pi 
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Figure 1: 25% (3 NSTX equilibrium. Adjacent con- 
tours are equidistant in 1/B. The thin, gray line is 
the separatrix and the cross marks the magnetic axis. 
The poloidal- field coils, carrying approximately 200 
kA of current each, are easily identifiable to the right. 



At NSTX, co-injection of an 80 kV deuterium 
beam is used to heat the plasma (co- injection is 
favored over counter-injection because it results 
in lower prompt losses) ||. Note that the beam 
directions are relative to the plasma current, which 
is in the opposite direction of the toroidal magnetic 
field, and a co-injected ion thus has vii < 0. Here, 
however, we will first launch a counter-injected ion, 
born with its gyro center just below the midplane 
at Rn C = 0.70 m. The resulting orbit is shown in 
Fig. g. Plotting the time evolution of the normalized 
magnetic moment A [see Fig. [|, the adiabaticity is 
evident; A oscillates around the constant mean with 
the amplitude spiking when the ion passes through 
the minimum of B along the orbit. The longitudinal 
invariant J|| is practically constant. 

Next, we launch a co-injected ion with its gyro 
center at exactly the same location as before for 
the counter injection. As can be seen in Fig. || the 
solution to the equation of motion now takes on a 
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Figure 2: Quasi-periodic orbit of a counter-injected 
deuterium beam ion. The dashed line is the limitcr. 
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E = 80kV ( V || = 0.87v, R gc = 0.70m, Z gc = -0 . 1 m ) 



ml "in 



I! Illm, 



\ nil!' 



Q n / 27T X t 



III 

1 60 



Figure 3: Adiabatic time evolution of normalized magnetic moment A (thin, gray line), longitudinal invariant 
J|| (solid line), and gyro-averaged magnetic field strength (B) (dashed line) for the orbit shown in Fig. |[ 



completely different character, as the adiabaticity 
condition is violated. The variation of the magnetic 
field along the orbit is actually less than it was for 
the counter- injected ion [see Fig. |j; but because the 
gyro period is longer and the drift period is shorter, 
the change over one gyration becomes sufficiently 
large to destroy the adiabaticity of both A and Jh . 
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Figure 4: Chaotic orbit of a co-injected beam ion. 



We now move the gyro center five centimeters 
outward, to the edge of the magnetic well. In the 
magnetic well the B contours are nested, closed loops 
and roughly overlap the flux surfaces. The variation 
of B along a flux surface is thus greatly reduced. 
Moving closer to the magnetic well, the adiabaticity 
parameter will gradually decrease, and launching the 
ion at R gc = 0.75 m, the character of the motion 
is again changed [see Fig. [|. Notice how the ion 
shifts between three different quasi-periodic orbits. 
This behavior can be understood by inspecting the 
time evolution of A and Jn [see Fig. 0. As can be 
seen, A does not oscillate around a constant average 
and is thus non- adiabatic. Instead it jumps between 
three different gyro averages, explaining the peculiar 
ion orbit. In contrast, Ju oscillates with a period 
of three drift periods around a constant average 
and thus qualifies as an adiabatic invariant. This 
phenomenon, with non-adiabatic A but adiabatic J\\ 
is known as super-adiabaticity Note that in the 
super-adiabatic regime, the longitudinal invariant J|| 
replaces the magnetic moment as the third constant 
of motion, and the motion remains integrable. The 
confinement of super-adiabatic ions should thus be 
similar to that of adiabatic ions. As the adiabaticity 
parameter is further reduced and we leave the 
super-adiabatic regime, however, only two constants 
of motion will remain, and the confinement will be 
drastically decreased as the motion becomes chaotic. 
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Figure 5: Non-adiabatic time evolution of A and Jn for the orbit shown in Fig. |[ A exhibits characteristic, 
erratic jumps as (B) passes through local minima. 
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Figure 6: Almost quasi-periodic orbit of a super- 
adiabatic co-injected beam ion on the edge of the 
magnetic well. 

The non-adiabatic beam ion shown in Fig. |] 
intersects the limiter, and it would have done so even 
if it had been adiabatic. This seems to be typical 
for NSTX: the non-adiabaticity only affects 80 kV 
beam ions already following loss orbits and should 
thus not have any substantial, negative effect on the 
confinement. The majority of the confined beam ions 
are super-adiabatic, similar to the orbit in Fig. |^. In 
fact, even orbits deep inside the well have a slightly 
non-adiabatic magnetic moment. 



In START, co-injection of a 30 kV hydrogen beam 
was found to have ~60% absorption efficiency p"2| . 
Based on numerical simulations, these beam ions 
were claimed to be strongly non-adiabatic. This 
would seem inconsistent with the good confinement 
and we note that Fig. 3a of Ref. [jjjj, which shows 
the time evolution of the magnetic moment for a 
typical beam ion and is said to demonstrate the 
non-adiabaticity, in fact indicates super-adiabatic, 
or possibly even fully adiabatic, beam ions. 

If the beam energy is increased to 110 kV in 
NSTX, the situation remains essentially the same: 
the vast majority of the beam ions that are not 
promptly lost, are super-adiabatic, and their mo- 
tion is consequently bounded. However, in future, 
burning-plasma experiments, such as the VNS p0[ , 
the fusion-born a-particles will display the whole 
range of dynamics: adiabatic, super-adiabatic, and 
non-adiabatic. As an example, we show the time 
evolution of the magnetic moment for an alpha, 
born just below the midplane and three centimeters 
outside of the magnetic axis, for a 38% (3 VNS 
equilibrium [see Fig. [g]. In addition to possible direct 
non-adiabatic losses when a giant jump in A causes 
transition to a loss orbit, there is also an indirect loss 
mechanism. As can be seen, the jumps in A force 
the a back and forth across the trapped-passing 
boundary. When driven by ion-ion collisions, this 
is of course the process that generates neo-classical 
transport. For fusion-born a-particles, the drift 
period is many orders of magnitude shorter than the 
collisional timescale, and the non-adiabatic transport 




Figure 7: Non-adiabatic time evolution of A but adiabatic time evolution of J|| for the orbit shown in Fig. |6| 
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Figure 8: Non-adiabatic time evolution of A and J|| for fusion a in VNS. 

mechanism described above can thus be expected to 
affect the overall a confinement. Yet another effect 
to consider is the Arnold diffusion Jl3| , which might 
be important if the magnetic-field ripple is large 
enough to destroy the axi-symmetry. 
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